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Abstract 

The Feynman all-coupling variational approach for the polaron is re-formulated and extended using the Hamiltonian 
formalism with time-ordered operator calculus. Special attention is devoted to the excited polaron states. The energy 
levels and the inverse lifetimes of the excited polaron states are, for the first time, explicitly derived within this all- 
coupling approach. Remarkable agreement of the obtained transition energies with the peak positions of the polaron 
optical conductivity calculated using diagrammatic quantum Monte Carlo is obtained. 
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Several studies have been devoted to the search of 
a Hamiltonian formalism equivalent to Feynman's path 
integral approximation to polaron theory. Bogolubov 
reproduced the Feynman result for the polaron free 
energy Isj] using time-ordering T-products . Yamazaki 
101 introduced two kinds of auxiliary vector fields to de- 
rive Feynman's ground state polaron energy expression 
with the operator technique, however he found no proof 
of the variational nature of this result. Cataudella et 
al. Istl formally re-obtained Feynman's polaron ground- 
state energy expression by introducing additional de- 
grees of freedom, but again their result could not be 
proved to constitute an upper bound for the polaron 
ground state energy. 

The study of the excited polaron states is of interest 
i. a. for its application to the polaron response prop- 
erties. In iQl a path-integral based response-formalism 
was introduced that was applied to derive polaron opti- 
cal absorption spectra in ||7[]. The results for the polaron 
response obtained in |01 were re-derived with a Hamil- 
tonian technique (Mori- formalism) in 1^. 

To the best of our knowledge, no explicit description 
of the polaron excited states has been derived within the 
"all coupling-" Feynman approach. Only for the limit- 
ing cases of weak and strong coupling approximations 
(and for a ID-model system) such excitation spectra 
were derived iM, llOll . 
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In principle, the spectrum of the polaron excited 
states can be derived indirectly - using a Laplace trans- 
form of the finite-temperature partition function. How- 
ever, it is not clear how to realize this program in prac- 
tice. 

The polaron excitation spectrum is interesting by it- 
self. E.g. the existence and the nature of "relaxed 
excited states", "Franck-Condon states", "scattering 
states" is understood from the mathematical structure 
of corresponding eigenstates. 

In the present letter we first present a re-derivation 
of the original Feynman variational path integral po- 
laron model [3] for the ground state, using a Hamil- 
tonian formalism, and we do provide a proof of the 
upper bound nature of the obtained ground state en- 
ergy. Furthermore, using Feynman's (Hamiltonian-) 
time-ordered operator calculus (and an ad hoc unitary 
transformation) we obtain explicitly - and for the first 
time - the excited polaron states that correspond to the 
Feynman polaron model. 

The novelty of the present approach consists (a) in 
the direct calculation of the energies and the lifetimes 
of the excited polaron states (within a Hamiltonian all- 
coupling approach - developed in this work - equiva- 
lent to the Feynman path integral polaron model) and 
(b) in the extension of the Feynman variational tech- 
nique to non-parabolic trial potentials. Although the 
time-ordered operator calculus is formally equivalent to 
the path-integral formalism, it is not obvious how to di- 
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rectly calculate the excited polaron states using path in- 
tegrals. 

The present work, formulated with the (Hamiltonian) 
time ordered operator calculus, thus provides an (equiv- 
alent) tool complementary with respect to the Feynman 
path integral approach to the polaron, to study the po- 
laron problem. Additionally we directly study the ex- 
cited polaron states. 

Consider an electron-phonon system with the 
Frohlich Hamiltonian 
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Here, the Feynman units are used: h- I, the band mass 
nib - I, the LO-phonon frequency wlo = 1 ■ 

The polaron partition function after exact averaging 
over phonon states is 
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where p = -j^. The "influence phase" of the phonons 
O [r^] in the the time-ordered operator calculus has the 
same form as in the path-integral representation. The 
polaron free energy is determined as 
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The trial Hamiltonian describes the electron interact- 
ing with a fictitious particle of the mass nif through an 
attractive potential V/: 
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The trial potential Vf is, in general, non-parabolic. The 
parabolic potential with frequency parameter w corre- 
sponds to the Feynman polaron model. 

Consider the "extended" partition function of the 
electron-phonon system 

Ze.vf - ZfZpol (7) 

where Zf is the partition function of a fictitious particle. 
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with Hamiltonian 
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The polaron free energy is expressed as the diff'erence 

Fpoi^F,„-Ff, (10) 

where is the free energy of the fictitious particle con- 
fined to the potential ^(r/). The free energies F^xt and 
Ff are determined similarly to (|5]l, with corresponding 
partition functions. In the zero-temperature limit, the 
free energies Fpoi, Fext and become, respectively, the 
ground-state energies E^^^_^j, E^^„ and E^^. 

The key element of the present approach is the unitary 
transformation 

t/ = e"'P/'. (11) 

Application of this canonical transformation results 
in the transformed "extended" Hamiltonian H'^^, - 
UHex,U-\ 
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This Hamiltonian can be represented as a sum of an un- 
perturbed Hamiltonian 



Ho = H,r + Hph 

and an interaction term 
1 
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Further we use the variational principle for the 
ground-state energy in terms of the time-ordered oper- 
ators following Ref. fll ill . The exact ground state |0) 
of the system with the Hamiltonian ( fT2b can be written 
in the interaction representation starting from the unper- 
turbed ground state |-oo): 

|0> = 'Z/(oo,-oo)|-oo) (15) 

where tlioo, -oo) is the time-evolution operator. 
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Here, (5 — » +0 and T denotes time ordering. 
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In the exact expectation value for the ground state en- 
ergy E'^^i = (q 0^, the phonons are eHminated us- 
ing the time ordered-operator calculus as in Ref. jl ill . 
The average of the interaction term becomes then 
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This means that the polaron ground state energy is ex- 
actly described using a retarded potential in the interac- 
tion representation, cf. Eq. (2.16) of Ref. ifTlll . 

The ground state energy satisfies the Ritz variational 
principle with a trial state. Choosing the trial state as 
the ground state of the Hamiltonian (fTST l. the variational 
principle can be written as 111 ill 



+ (cx, \T {1/,, (oo, -ex,) [H^^, (0) - Ho (0)]}| - «,) , 

(18) 

where 'Uir {oo, -oo) is the time-evolution operator corre- 
sponding to the trial Hamiltonian 

The exact polaron ground state energy is denoted here 
as (k), where k is the polaron translation momen- 
tum. We find an upper bound for E^ (k) substituting ([TTl i 
in (fTSl l and using the exact wave functions and energy 
levels of the trial Hamiltonian. The trial Hamiltonian 
(|6]l can be rewritten in terms of the coordinates (R,p) 
and momenta (P, n) of the center-of-mass and relative 
(internal) motions of the trial system with the masses 
M - \ + nif and ji - inf/(l + nifj using the frequency 
V = wM. The energy spectrum of the trial system is the 
sum of the translation- and oscillation contributions. 
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The eigenfunctions of the Hamiltonian (jSjl are products 
of translational- and oscillatory wave functions: 



l/'k;/,«,m (R, P) = -^e'^^'^ipin^yn (p) , (20) 



where </3/,„,,„ (p) is the 3D harmonic-oscillator wave 
function with a given angular momentum. The result 
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where v > w are the Feynman variational frequencies. 
The functional (l22l i can be reduced to the known Feyn- 
man result for the polaron ground-state energy. In the 
r.h.s. of (i22lat the polaron momentum k = 0, we intro- 
duce the integral over the Euclidean time: 



. = e ^ / dr. (23) 

+ vn' + 1 -'0 



After this, the summations and integrations in ( |25] ) are 
performed analytically, and we arrive at the Feynman 
variational expression for the polaron ground-state en- 
ergy: 
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The electron-phonon contribution in ( l22l i is struc- 
turally similar to the second-order perturbation cor- 
rection to the polaron ground-state energy due to the 
electron-phonon interaction (using states of the Feyn- 
man model iffk-jjijii as the zero-order approximation). 
Therefore we can estimate the energies of the excited 
polaron states when averaging the difference between 
exact and unperturbed Hamiltonians in (fTSl l with an ex- 
cited trial state. We then arrive at the following exten- 
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sion for the r.h.s. of ( l22b : 
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In the same approach, we obtain the inverse lifetimes 
for the excited states of the polaron: 
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The broadening of the excited polaron "non-scattering" 
states must be taken into account for an analytical study 
of the polaron optical conductivity. 

Using the above expressions, we determine the tran- 
sition energies for the transitions between the ground 
and the first excited state ?iQo^ie.i-e = £p - Ef\ Let 
us first consider the transition energies in which E'*/'^^'^' 
are calculated using optimal values of the parameters 
of the Feynman model obtained from the minimiza- 
tion of the variational ground-state energy Ef \ This 
method formally leads to the Franck-Condon (FC) ex- 
cited states, with the "frozen" phonon configuration cor- 
responding to the ground state of the polaron. Note 
that the existence of Franck-Condon states as eigen- 
states of the Frohlich polaron Hamiltonian has not been 
proved: Ref [10] suggests their non-existence as eigen- 
states for a simplified polaron model. Nevertheless the 
Franck-Condon concept can be significant, e. g. for ap- 
proximate treatments using a basis of Franck-Condon 
states, as indicative for the frequency of the maxima of 
phonon-sidebands, etc. 

In Fig. 1, the FC transition energies calculated with 
the approach introduced in the present work for polaron 
momentum k = are plotted as a function of the cou- 
pling constant a. They are compared with the peak 
energies of the polaron optical conductivity calculated 
using the diagrammatic Monte Carlo method (DQMC) 
lll2[ll3ll and with the peak energies attributed to polaron 
"relaxed excited states" (RES) in Ref. 101 ("DSG"). 
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Figure 1 : Franck-Condon transition energies as a function of the cou- 
pUng constant compared to the lowest-energy peak position of the po- 
laron optical conductivity from Ref. 01 and the maximum of the po- 
laron optical conductivity band from Ref. fT3l . 



The DQMC and DSG main-peak energies are close to 
each other in the whole range of the coupling strength. 
In the range 4 ^ a ^ 10, the present result for the 
transition energy is close to the DQMC and the DSG 
peak energies. Furthermore, in this range of a, the non- 
monotonous behavior of the curvature is remarkably the 
same for the DQMC and DSG peak energies and for the 
present result. 

There is a remarkable agreement between the peaks 
attributed to the RES in Ref. 101, the peak positions 
obtained within the strong-coupling approach, Eq. (3) 
of Ref. 11311 . and the positions of the maximum of the 
optical conductivity band calculated in Ref. lll2ll using 
DQMC. It is reasonable that the three aforesaid peaks 
must be interpreted in one and the same way. In order 
to clarify this, we can refer to Ref. 11211 . In the strong- 
coupling regime, the dominant broad peak of the po- 
laron optical conductivity spectrum can be considered 
as a "Franck-Condon sideband" of the "groundstate to 
RES-transition", even if this latter transition can have a 
negligible oscillator strength (see also |[141 ). The optical 
conductivity spectra of Ref. 11311 in the strong-coupling 
approximation have been calculated taking into account 
the polaronic shift of the energy levels. The pola- 
ronic shift in Ref. ifTill has been calculated with the 
Franck-Condon wave functions (i. e., with the strong- 
coupling wave functions corresponding to the "frozen" 
lattice configuration for the ground state). Note that the 
exact excitation spectrum of the Frohlich-Hamiltonian 
might be devoid of Franck-Condon eigenstates, cf. Ref. 
ildH ). It should be remarked that the maxima of the FC- 
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sideband structures of Ref. |01 are positioned at the fre- 
quency Q - v, i. e., at the transition frequency for the 
model system without the polaron shift. 

The Franck-Condon peak energies calculated in the 
present work also take into account the polaron shift. As 
follows from the above analysis, in the strong-coupling 
limit they must correspond to the Franck-Condon peak 
energies of the strong-coupling expansion of Ref. lll3ll . 
The agreement of the position of the maxima of these 
peaks with those attributed to transitions to the RES 
in Ref. 10] shows that in the strong-coupling range of 
a, the latter should be associated to the Franck-Condon 
sidebands rather than to the RES. 

Another approach, in which the parameters of the 
first excited state are determined self-consistently (Ref. 
1 1411). was used i. a. to calculate (in the strong-coupling 
case) the (lowest) energy level of the relaxed excited 
state (RES). The transitions from the polaron ground 
state to the RES correspond to a zero-phonon peak in 
the optical conductivity. 

For the study of the energies of excited states of the 
polaron, a variational approach requires special care, 
because the excited states of the polaron are not stable. 
A variational approach, strictly speaking, is only valid 
for excited states when the variational wave function of 
the excited state is orthogonal to the exact ground-state 
wave function. 

For the estimation of the energy of the first RES with 
our present formalism, we determine a minimum of the 
expression dZST l in a physically reasonable range of the 
variational parameters. In order to determine that range, 
we refer to Ref. jlsll, where the energy of the polaron 
RES is calculated variationally within the Green's func- 
tion formaUsm. 

The expression for the RES energy in Ref. ifisll con- 
tains the electron-phonon contribution corresponding to 
the second-order perturbation formula. It differs, how- 
ever, from the weak-coupling second-order perturbation 
expression by the choice of the unperturbed states: in 
Ref. 11 1511 those are variational states rather than free- 
electron states. There exists some analogy between our 
approach and that of Ref. I>15il . The latter, however, 
does not take into account the translation invariance of 
the polaron problem. 

In Ref. 1151], the energy of the polaron RES is cal- 
culated variationally. The unperturbed wave function of 
the RES is chosen orthogonal (due to symmetry) to the 
unperturbed ground state wave function. In the present 
approach, this orthogonality is also exactly satisfied be- 
cause of symmetry. 

The expressions for the polaron RES energy of Ref. 
1 151] contain singularities, which occur when the ener- 



gies of the unperturbed ground state and that of the first 
excited states are in resonant with the LO-phonon en- 
ergy. These singularities are related to the instability of 
the excited polaron with respect to the emission of LO- 
phonons. Using the same reasoning as in Ref. 11511 we 
search for a local minimum of the polaron RES energy 
in the range where the confinement frequency v of the 
Feynman model satisfies the inequality \' > 1 . The in- 
stability of the excited polaron state is then avoided. 

The resulting numerical values of the transition en- 
ergy to the first RES as a function of a are shown in 
Fig. 2. They are compared with the numerical-DQMC 
pea k energies of the polaron optical conductivity band 
il2lll3ll . with the FC transition energies obtained in the 
present work, and with the leading term of the strong- 
coupling approximation for the RES transition energy 
from Ref. IM]. 
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Figure 2: The transition energy for the transition from the polaron 
ground state to the first RES(solid black cun'e) and to the first excited 
FC state (dashed red cun'e) as a function of a obtained in the present 
work, compared with the maximum of the polaron optical conduc- 
tivity band from numerical DQMC {black squares, Ref. (T^). The 
dashed-dot green cun>e: the strong-coupling result for this transition 
energy as given in Ref. QJ. 



For a < 2.5, there exists no minimum of E' 



{le.xc) 
P 



in the 



range v > 1. We can interpret this result as a manifes- 
tation of the fact that for decreasing coupling strength, 
the RES is suppressed at sufficiently weak coupling. We 
see that for sufficiently small a (a < 6), the RES tran- 
sition energies show good agreement with the DQMC 
peak energies, what confirms the concept of RES devel- 
oped in Refs. ifTi ffl . For higher coupling strengths, the 
DQMC data appear to be closer to the FC (rather than 
to RES) transition energies. This result can be an indi- 
cation of the fact that with increasing a, the mechanism 
of the polaron optical absorption changes its nature as 
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suggested in Ref. II13I1 . from a regime with dynamic lat- 
tice relaxation (for which the RES are relevant) at weak 
and intermediate coupling to the Franck-Condon ("LO- 
phonon sidebands"-) regime at strong coupling. 

In summary, we have re-formulated the Feynman all- 
coupling path integral method for the polaron problem 
within a Hamiltonian formalism using time-ordered op- 
erator calculus. This reformulation allows us to de- 
scribe not only the free energy and the ground state, but 
also to directly determine - for the first time - the ex- 
cited polaron states that correspond to the Feynman all- 
coupling polaron model. A variational procedure for the 
polaron RES energy has been developed, within the for- 
malism presented in this work, which provides results 
i.a. in agreement with the strong-coupling limit of Ref. 



1 14|] . The present treatment offers the prospect of fur- 
ther elucidation of the nature of the polaron resonances 
("relaxed excited states" versus "Franck-Condon side- 
bands" [9]) at intermediate coupling. 
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